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Abstract 

We study the phase diagram of the three-dimensional SU(3)+adjoint Higgs theory 
with lattice Monte Carlo simulations. A critical line consisting of a first order line, 
a tricritical point and a second order line, divides the phase diagram into two parts 
distinguished by (Tr Afy = and 7^ 0. The location and the type of the critical line 
are determined by measuring the condensates (Tr Aq) and (TtAq), and the masses of 
scalar and vector excitations. Although in principle there can be different types of 
broken phases, corresponding perturbatively to unbroken SU(2)xU(l) or U(l)xU(l) 
symmetries, we find that dynamically only the broken phase with SU(2) xU(l)-like 
properties is realized. The relation of the phase diagram to 4d finite temperature QCD 
is discussed. 
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1 Introduction 



The motivation for studying the three-dimensional (3d) SU(3) + adjoint Higgs gauge 
theory is twofold. First of all, this case is interesting since the 3d SU(3) + adjoint 
Higgs theory is an effective theory of finite temperature QCD in the weak coupling 
domain |J. The requirement of small couplings means that this effective theory is 
accurate only in the limit T ^> Aj^g, not in the phase transition region T « T c 
For pure gauge SU(3) theory, this is related to the fact that the phase transition at 
T = T c has to do with the breaking of the Z(3) symmetry. This symmetry is lost in 
the effective theory, but some traces of it remain, as will be discussed below. 

The second interesting aspect of the 3d SU(3) theory is that it already has several 
of the properties of the SU(5) case, which is relevant for some GUTs at finite tem- 
peratures ||. These properties are not shared by the special SU(2) case, where the 
structure constants d abc vanish. Such properties are the existence, in perturbation the- 
ory, of different types of broken phases with the associated various gauge groups, and 
of the corresponding monopoles. However, the spectrum of various broken phases is 
of course richer in SU(5) than in SU(3), and the SU(5) action has also two non-trivial 
scalar self-couplings in contrast to the SU(3) action, which only has one. Nevertheless, 
we expect a similar qualitative behaviour. 

Some conjectures concerning the phase diagram of the 3d SU(3) + adjoint Higgs 
theory were put forward in 0. The purpose of this paper is to determine the phase 
diagram numerically. 

The actions of 3d SU(iV < 3) + Higgs theories normally depend on two variables^]: 
the scalar mass m\ and the scalar self-coupling A3. On the mean field level the system 
has two phases: a "symmetric" phase at m\ > and a "broken" phase at m\ < 
at any A3. The problem now is to determine the phase diagram in the full quantum 
theory. This question has previously been answered with numerical lattice Monte Carlo 
simulations for a number of theories: SU(2) + fundamental representation Higgs H- 
0, SU(2)xU(l) + fundamental Higgs |2j, and SU(2) + adjoint Higgs p], ||. The 



case of U(l) + fundamental scalar Higgs (the Ginzburg-Landau theory) has also been 
extensively studied JTJ| . 



For small Higgs self-coupling, A3 <C g%, the tree- level second order transition is, in all 
cases studied, radiatively changed into a first order transition. Its strength decreases 
with increasing A3. The central and often very difficult question is what happens at 
larger A3: does the first order line terminate (so that the two phases are analytically 
connected) or is there a phase transition line across the whole phase diagram (so that 
there is an order parameter vanishing in one of the phases). In this respect, the phase 
diagrams for the different symmetry groups differ in quite interesting ways from each 
other. A good illustration of the difficulty of the problem is the fact that the nature 
of the phase diagram in the superficially simplest case, the Ginzburg-Landau theory, 



D As discussed below, in some cases there can in principle be more couplings than just two. 
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is not yet conclusively determined in the large A3, extreme type II, regime. 

For SU(3) + adjoint Higgs theory we shall see that the phase diagram is divided into 
two parts by a phase transition line which contains a first order line, a tricritical point 
and a second order line. In contrast, the SU(2) + adjoint Higgs theory is observed to 
have no transition at large A3 |13|]. The reason why it is possible to have a qualitatively 
different behaviour in SU(2) and SU(3) is that for all SU(iV > 3) there is a gauge- 
invariant local order parameter (Ty Aq) sensitive to the breaking of the Aq <-» —Aq 
symmetry of the theory. 

The plan of the paper is the following. In Sec. |2] we define the theory in contin- 
uum and on the lattice. In Sec. § we present some perturbative estimates for the 
different observables measured, and in Sec. ^ we briefly review the relation to 4d finite 
temperature QCD. The simulation results are in Sec. §, and the conclusions in Sec. ||. 



2 Definition of SU(3) + adjoint Higgs theory 

The theory we study is defined by the following super-renormalizable Lagrangian: 



L[AIAI\ = ^F^ + TT[D t ,A ][D l ,A ]+mlTTAl + \ A (TTA 2 ) 2 , (1) 

in standard notation (Di = di + ig^Ai). The notation Aq = AqT 11 for the adjoint scalar 
is chosen with its origin from dimensional reduction of QCD in mind. For SU(3) and 
SU(2) Tr^4g = |(Tr Aq) 2 so that only one scalar self-coupling appears. In principle, 
there could also be two additional super-renormalizable couplings h 3 , h 5 appearing as 

6L = haTrA% + h s TrA%TrA%. (2) 

However, these terms do not arise in the dimensional reduction of finite temperature 
QCD, and thus we assume that h 3 = h 5 = 0. 

In the absence of h^,^, the theory in Eq. flip is symmetric under A — > —Aq. 
This symmetry was called "R-parity" in [|7]]. In terms of 4d physics, this symmetry 
is related to the usual discrete transformations CT, P |I5|, [TjJ. However, it should be 



clearly stated that the breaking of the Aq — ► — A Q symmetry to be discussed below, 
does certainly not imply spontaneous breaking of any of the discrete symmetries of 
finite temperature QCD, since the broken phases are not physical from the point of 
view of QCD @. 

Due to super-renormalizability, of the couplings in Eq. ([I]) only the scalar mass m\ 
is scale dependent. In the MS scheme 



m ^ ) = /2 = 20(3^ 2 - A A )A A , (3) 
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where A is a constant of dimension GeV. The dynamics of the theory is thus determined 
by the dimensionless variables 



i/3 i/3 

where g§ has been used as a natural mass unit. 

Instead of regulating the theory using the MS scheme one can as well use the lattice 
scheme with a lattice constant a. The two schemes have to be matched so that they, in 
the limit a — > and for given gf, y, x, give the same physical results. Due to the super- 
renormalizability this only requires tuning the bare mass term of the lattice action. 
Introducing the gauge coupling constant Pq by 

Pg = A> (5) 
the lattice action S = S[Ui(x), A Q (x)} becomes 

x i<j 6 

12 

- — £ Tr A (x) U t (x) A (x + i) u} (x) 
+ 2 |/3 2 Tr^(x) +x^[Tr^(x)] 2 | , (6) 
where -Pij(x) is the plaquette and where the coefficient of the quadratic term is [17] 



36 f 6 ,„ _ ,3.1759115 
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A = 1 + ^-(6 + 10 ^^ r 



, [(60a;- 20x 2 )(ln/3 G + 0.08849) + 34.768a; + 36.130] L (7) 

Remarkably, the matching of lattice and continuum can be carried out analytically even 
including terms of order a [18|. This implies that the gf, y, x in Eq. (H) are modified by 
corrections of order 1/Pg- The additive corrections for y have not yet been computed: 
the 0(a) correction is, for dimensional reasons, ~ ag®f(x,y) and a 3-loop computation 
would be needed. For the two other parameters the improvements are 

2 / 1.994833\ 2 

[93 j improved — I 1 H T. I 9zi 



G 



^improved = x - -J-(0. 328432 - 0.835282a; + 1.167759a; 2 ), 

Pg 
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where CJ^jX £1X6 the parameters appearing in Eqs. (H) - ©- 

The improved relations between the condensates of the scalar field Aq in the lattice 
action in Eq. (||) and in the continuum areQ 



<Tr^> 



cont 



3 



3.1759115^ G 



2tt 



6lT 

2 (ln/5 G + 0.6678), 



3.409891 - 0.729850x\ 
Jg J 



(TrA*) 



cont 



9l 



<Tr^> 1- 



5.114364 



13, 



G 



0.437910x\ 



(9) 
(10) 



In Eq. (Q) an 0(ag^) additive 3-loop correction is still not known. The l//?G-terms are 
numerically quite large at the value (3g = 12 we have used in practice (up to ~ 40%), 
and implementing them brings the lattice results significantly closer to the perturbative 
results (deep in the perturbative regime). 



The MS-scheme regularized operator (Tr A^) cont is, in fact, scale dependent |B|, and 
has been defined at the scale \i — g\ in Eq. (|9|). The scale dependence arises at the 
2-loop level and is, for general N, 



(TrAW) = (Tr^) + 



In terms of the effective potential (note that (TrAQ) cont = dV (mm) / dm^) , the scale 
dependence arises from the graph (solid line = Aq, wavy line = A{) 




X ^-1)N / _g_ 3 
16tt 2 3 V 2m 3 4 



(12) 



It is amusing to note that this term is precisely the 0(g hi(g)) Toimela term E0, 131 in 



the free energy of the QCD plasma. 

The gauge invariant operators of lowest dimensionalities in the action defined by 
Eq. (^|) are as follows: 

• Dim=l: TrAg, 

• Dim=3/2: TrAg, 

• Dim=2: (Tr A 2 ) 2 , e ijk Tr A F jk , 

• Dim=5/2: e ijk Tr A 2 F jk , 

• Dim=3: Tr FijF k i, Tr [Di, Aq] [Dj, Aq] (with various spin projections). 

These operators can be used for mass measurements in the different quantum number 
channels. 



7 We thank Guy Moore for providing us with the improved expression of (Aq). 
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3 Perturbation theory 



Because of confinement, perturbation theory is not well convergent in the symmetric 
phase of the theory and is thus, in general, of limited usefulness in the study of the 
phase structure. Nevertheless, perturbation theory does work in the limit of small 
x = Aa/(?3 when the transition becomes very strong, and it is worthwhile to go through 
its predictions there. 

Let us first fix a gauge and parametrise a constant diagonal SU(3) background field 
as follows: 



5,T 3 



B»T S 



( q + p 
q — p 
-2q 



with 



Tr B 2 = 6q 2 + 2p 2 , Tr B 3 = 6q(p 2 - q 2 ) 



Tr£ 4 



1 



(Tr 5 



2\2 



(13) 



(14) 



(Note that Tr B 2 = g 3 2 Tr (A ) 2 ^ g 3 2 (TrA 2 ).) If p = or p = ±3q, then the SU(3) 
symmetry is broken to SU(2) xU(l), otherwise it is broken to U(l)xU(l). These are 
perturbative statements; the only symmetry that can be broken in the full quantum 
theory is A «-> —A , which is signalled by (TtAq). 



3.1 The x -> limit 

Let us compute the 1-loop effective potential V\(x, y; q, p) in the background in Eq. fll3"D. 
One finds, in the Landau gauge, 

2\2 



g 3 b V(x, y; q, p) = y(6q 2 + 2p 2 ) + x(Qq 2 + 2p 2 ) 

-7^(|3g + p| 3 + |3g-p| 3 + 8b| 3 ) 
Sir 

\-{[y + 6x(6g 2 + 2p 2 )f 2 + 7[y + 2x(6q 2 + 2p 2 )f' 2 } . 



12n 



(15) 



The scalar loop terms (last line in Eq. fjl5|) ) become negligible when x — > (apart from 
a constant) and we shall neglect them to begin with. Then, for p — 0, 

1 ^ 2 



^3 V(x,y;q,0) = 3Qxq 



Attx 



+ 



6.r 



y - 



8tt 2 x 



Thus, for 



y = y, 



Hoop 



(X) 



Itt 2 x 



the system has two coexisting states (a first order transition) with 

1 



Asymmetric > Abroken 



Atix 



(16) 



(17) 



18) 



One stable and one metastable state exists for 



y _(x) = < y < y+(x) 
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■Vc{x). 



(19) 



These results for p = can now be extended to the whole (g, p)-plane. Indeed, 



V 1 (x,y;q> 0,p = 3q) = V x (x, y; 2q, 0), Vi(x,y;q< 0,p = -3q) = V(x,y;2q,0). 



This implies that, at y — y c {x), the system has the following degenerate potential 
minima (Fig.|I]): 

(0) a symmetric minimum at q = 0, p = 0, 

(1) a broken minimum at q = l/(87rx), p = 3/(87nc), 

(2) a broken minimum at q = l/(4irx), p = 0. 

Moreover, in the parametrisation in Eq. flT3| ) one has the freedom of permuting the 
diagonal elements and the potential has to be invariant under the transformations of 
q,p corresponding to these permutations. One sees that the fundamental region, which 
determines the potential over the whole plane, can be chosen to be that bounded by 
the two lines p = 0, p = 3q. Thus, for each broken minimum there are two more 
minima corresponding to cyclic permutations of the diagonal elements. All these min- 
ima correspond to breaking to SU(2) xU(l); there are no local minima corresponding 
to breaking to U(l)xU(l) (a U(l)xU(l) minimum would require that all the vector 
masses cubed appearing on the 2nd line in Eq. (|T5| ) are non-vanishing, but this is not 
the case in any of the minima considered). 

To get rid of the unphysical extra symmetries related to the permutation of the 
diagonal elements of B in Eq. (|T3|), it is useful to write the 1-loop potential in terms of 
Ty B 2 and Tr5 3 . Effectively, one inverts the cubic Eqs. (|14D and inserts into Eq. (p!5|). 
The result is 



This form, containing the solution of a cubic equation, is not very transparent, but it 
actually is quite simple, as shown in Fig.|l]. Now only the genuinely different minima, 
one symmetric and two broken ones, appear. 

Summarising, for x — > perturbation theory predicts that the system has the fol- 
lowing three phases: 



(20) 



Vt{x, y; Tr B 2 , Tr B 3 ) = yTr B 2 + x(Ti B 2 ) 2 




(21) 
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TrB 2 A 



•© 




8nx 




TrB 2 = [6(TrB3)2] 1/3 



3 TrB 3 
327t3 x 3 



Figure 1: The minima of the 1-vectorloop potential on the critical line y = y c (x) 
(Eq. (tL"7j) ). Left: The minima in the (g,p)-plane. Those with the same number corre- 
spond to permutations of the eigenvalues of B in Eq. (p~3|) . The fundamental region is 
bounded by the two lines p = 0, p = 3q. Right: The minima in the (Tr B 2 , Tr i? 3 )-plane. 
The boundary curve is the map of the triangle in the left panel. 



• one symmetric phase with TrB 2 = TrB 3 = for y > yl loop (x), 

• two broken phases for y < y\ p (x), distinguished by the sign of Tr B 3 . 
In the broken phase at the critical point y = y c , 



TrB 2 
TrB 3 



iir 2 x 2 



3.80 



0.1 



x 



32lT 6 X 6 \ X 



(22) 



3.2 Condensates 

Using the effective potential in Eq. (fLSD , one can calculate the 1-loop perturbative 
approximation for (TrA 2 ), 



gf^rA 2 ) = g 3 ^ = 6q 2 - J- Uy + SGxq 2 ^ 2 + 7(y + I2xq 2 f' 2 
ay 87r L 



(23) 



where q is obtained by minimizing the potential in Eq. (|TJ) with p = and is zero in 
the symmetric phase. A similar calculation yields 



/v'nr.i;;) : U»r 3,/ 



57T 



(y + 36xq 2 ) 1/2 - (y + 12xq 



2U/2 



(24) 



Figure 2: Schematic phase diagrams of systems with a tricritical point (left) and a 
quadruple point (right), shown in a 3-dimensional coupling space (x,y,h), where h is 
an external field which couples to Tr Aq. The surfaces are 1st order phase transitions, 
and dashed lines 2nd order (Ising) transitions. The 1st order surfaces join along lines of 
triple points, where three phases are separated by 1st order phase transitions. At the 
quadruple point (right) four phases are separated by 1st order transitions; this is an 
intersection point of four lines of triple points and six surfaces of 1st order transitions. 



At the limit x -> considered in Eq. fl22|), g^ 2 (TiA 2 ) = Tr B 2 - (y 1/2 /ir), ^ 3 (Tr A 3 ) = 
TyB 3 . Note that the perturbative MS value of (TtAq) is negative in the symmetric 
phase. 



3.3 What happens at larger xl 

When x increases, the scalar contributions become more important, perturbation the- 
ory becomes less accurate, the three well separated phases in Fig.|l| approach each other, 
and the transition gets weaker. In the case of the SU(2) + adjoint Higgs theory this 
leads to an endpoint of the first order line fpT3|| : there is no local order parameter, and 



the phases are believed to be analytically connected. In the present case, the transi- 
tion also becomes weaker with increasing x. However, now there is a gauge invariant 
local order parameter, (TrA^), which can signal the breaking of the Z(2) symmetry 
A — > —A of the theory. Thus we expect that the (x, y)-plane phase diagram can be 
disconnected by a critical line containing a first order transition at x < x c , a tricritical 
point at x = x c , and a second order line at x > x c . 

However, as is standard in the case of tricritical transitions, in order to see the full 
phase diagram one has to use three couplings/fields. In our case the couplings are 
x, y and h, an external field with a coupling hTr Aq (see Eq. ([|)). The schematic 
3-dimensional phase diagram is shown in the left panel of Fig. |^. We note that in this 
full coupling space the three phases are still analytically connected through the large 
x region. 
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In principle, more complicated phase structures are also possible. At large enough x 
a new phase with (Tr Aq) = 0, (Tr Aq) > may appear (compare with the symmetry 
breaking patterns shown in Eq. (|14]) !). There is no local order parameter separating 
this phase form the symmetric one, but perturbatively it corresponds to a phase with 
the gauge group U(l)xU(l). If this occurs, there will be a quadruple point at some 
value of x, where four phases separated by 1st order transitions can exist (see the right 
panel of Fig. ||). The phase structure suggested in belongs to this class, even though 
the authors discussed the phases only in the (x, y)-plane. 

This kind of phase structure with a quadruple point is known to exist in some 
theories, for example, in the 3d 3-state Potts model and, indeed, in the finite 
temperature pure SU(3) gauge theory in 3+1 dimensions. These models have an exact 
Z(3) symmetry, and they have 3 degenerate broken phases and 1 unbroken phase, which 
can all exist at the quadruple point. For the case of SU(3) gauge theory, the 3d phase 
diagram is spanned by 1/T and external fields coupling to the real and imaginary parts 
of the Polyakov line. Indeed, the phase space of the Polyakov line in the complex plane 
does resemble the phase structure of the SU(3) + adjoint Higgs theory, shown in Fig. |l|. 
If one considers the SU(3) + adjoint Higgs theory as a dimensionally reduced version 
of the 4d gauge theory, it is appealing to think that the phase structures of the two 
theories could be similar. However, one has to bear in mind that in the 3d adjoint 
Higgs theory the 3-fold symmetry of the original 4d theory is strongly broken, as will 
be discussed in the next section. Even then it is not impossible that the Z(3) symmetry 
is dynamically generated in the close proximity of the would-be tricritical point (this 
can be compared with the dynamically generated Ising symmetry in the SU(2) gauge 
+ fundamental Higgs theory |TT[). If this truly happens, the would-be tricritical point 
could transform into a quadruple point. 

The question about the nature of the phase diagram is settled numerically below, 
and it will turn out that the phase diagram is the standard tricritical one, as shown in 
the left part of Fig. |2|. 



4 Relation to 4d and to Z(iV) symmetry 

The relation between x, y and 4d physics (T, Aj^g, N, Nf) is worked out explicitly to 
leading + next-to-leading order in [Q, [22] . For N = 3, Nf = the answer is 

x = ^ ^ 2(47re _ 7£ _ 3 /ll r) 

8tt 2 

1 , (25) 



11 ln(5.371T/A 



MSJ 



V = y d r(x) = -^ + -^ + 0(x). (26) 

87l z X l07T z 

As discussed, the 4d finite T theory without matter in the fundamental representation 
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x Volumes 



0.10 
0.15 
0.20 
0.25 
0.30 



12 3 

16 3 

16 3 , 

16 3 , 

16 3 , 



32 3 

24 3 , 32 3 , 48 3 
32 3 



Table 1: The lattice sizes used at (5g = 12. For 0.15 < x < 0.30 several values of y were 
used for each x, both above and below the transition. For x < 0.20, multicanonical 
simulations were used at y = y c (x). 

possesses an extra symmetry, the Z(iV) symmetry, for which the thermal Wilson line 



is an order parameter: if Zk = exp(i2nk/N) one can define a transformation of the fields 
so that the action is invariant but L — > z^L. Thus (L(x)) acts as an order parameter. 

This symmetry, however, is lost in the reduction process. The reason is simply that 
in the reduction process mass scales ~ 7rT are integrated over. However, under a 
Z(N) transformation a field configuration with small fields is transformed to one with 
A ~ nT/g and the reduction process cannot accurately represent such large scales. 

Even though the exact symmetry is lost, part of it is restored by radiative effects 
in the effective theory. The leading order term in Eq. (|26|) coincides with the leading 
order critical line Eq. (|TTD. This can be seen as a remnant of the Z(N) symmetry, since 
in this approximation, the effect of the Z(N) transformations of the 4d theory is to 
move the system from one of the three degenerate minima of Fig. [T] to another. This 
degeneracy is lost when higher-order corrections are taken into account, and the true 
transition line of the 3d theory does not agree with the dimensional reduction line. 
Moreover, in the effective 3d theory, only the vertices marked by 1 and 2 are related 
by the symmetry of the theory; in hot QCD (and in the 3d 3-state Potts model) all the 
vertices are equivalent. It is thus also natural that in the effective theory the middle 
point of the triangular region in Fig. [1] plays no special role, in contrast to the situation 
in hot QCD, where the middle point corresponds to the confined phase. 

5 Simulation results 

The primary aim of the simulations is to resolve the nature of the phase diagram and 
find the (x, j/)-plane critical curve y = y c (x). We also measure the screening masses 
(inverse screening lengths) on both sides of the transition line at various values of x. 




(27) 
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The simulations were performed using (5q — 6/(fl'f a ) — 12, with x in the interval 
0.10 < x < 0.30. The use of only one lattice spacing (one (3g) precludes the extrap- 
olation of the results to the continuum limit. However, we expect the finite lattice 
spacing effects to be small enough for the purpose of mapping the phase diagram. This 
is supported by our experiences from the measurement of the phase diagram of the 
SU(2) + adjoint Higgs model |§. 

The lattice volumes and the values of x used in this study are listed in Table [I]. For 
each x and lattice size, several runs at various values of y were performed. The total 
number of runs was 72, with approximately 2.1 node- years of cpu-time on a Cray T3E. 
This corresponds to ~ 2 x 10 16 floating point operations at 130 Mflops/node. 

5.1 Local order parameters 

For x < 0.2, the transition point y = y c (x) was determined with multicanonical sim- 
ulations using only lattices of size 12 3 and 16 3 . The transition here is so strongly of 
the first order - that is, the latent heat and the surface tension are so large - that the 
system tunnels from one phase to another too infrequently, even when multicanonical 
simulations are used. The value of y c {x) was determined to be the value of y where 
the probability weight of each of the 3 phases is equal. 

The strength of the transition is illustrated in Fig. ^, where we show the probability 
distribution in the (TrAjj, Tr A^-plane on the transition line for x = 0.15. ..0.25. For 
small x, the three peaks are very strongly separated, but when x approaches 0.25, the 
peaks join. It should be noted that in each case, the two Tr Aq > peaks are always 
connected by a "tunnelling channel" to the TrA^ pa phase. When x = 0.15, the 
relative probability density in these channels is suppressed by a factor ~ 10~ 9 when 
compared with the peaks, making them utterly invisible in Fig. |3], plotted with a linear 
scale. At x = 0.2, this suppression is "only" a factor of ~ 0.01, and at x = 0.25, there 
is no significant suppression any more. Note that the magnitude of the suppression is 
not universal and it depends very strongly on the volume of the system. Indeed, in 
large enough volumes the magnitude of the suppression can be related to the interface 
tension a between the phases. However, the volumes used here are too small for a 
reliable determination of a. 

It should be noted that there is no "tunnelling channel" directly connecting the two 
Tr^Q > peaks even near the tricritical point. Thus, the tunnelling from one of 
these peaks to the other always proceeds through the region of the Tr pa peak. 
This indicates that there is no dynamical generation of the Z(3) symmetry, and the 
behaviour is the standard tricritical one. 

In fact, we did not observe any sign of the quadruple point and the associated new 
phase in any point of the (x, y)-plane (since (TrAg), (Tr A^) are always simultaneously 
small or large, see below). This indicates that the phase diagram shown in the right 
panel of Fig. is not relevant here. 
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Figure 3: The probability distributions in the (Tr Aq, Tr A§)-plane at y — y c (x), x = 
0.15, 0.20, 0.25. The data is from 16 3 , (3g = 12 simulations; the simulations at x = 0.15 
and x = 0.20 are multicanonical. 

In Fig. |] we show the behaviour of the dim-1 and dim-(3/2) condensates (scaled by 
a proper power of g% to make them dimensionless) , 

<T^E T ^o/S3 2 > and <|^£Tr^o/S3 3 |>, (28) 

at 0.15 < x < 0.30 across the transition. Note that in the latter case a projection to 
positive values is needed, since otherwise the operator would always yield zero when 
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Figure 4: The condensates (-^3 Tr A^gf) and (\^ s J2 x Tr Al/ gfl) measured across 
the transition at x — 0.15, 0.20, 0.25 and 0.30 (converted to continuum units accord- 
ing to Eqs. (H), (|ToD ). The broken phase values compare quite well with the 1-loop 
perturbative ones in Eqs. ([£§) for x = 0.15. 



measured from a finite volume. In the infinite volume limit, J2x Tr ^o/dH) is 
perturbatively given by Eq. In the figures, the lattice operators in Eq. fl2"8|) are 



converted to continuum units according to Eqs. @, (]10f) . 

At x < 0.2 the transition is strongly of the 1st order, and there is a substantial meta- 
stability range across the transition point y c . However, the operator ( | -Ag J2 X Tr A^/g^l), 
which has no additive renormalization, remains a good order parameter at all values 
of x, even when the metastability disappears at x>0.25. Its deviation from zero 
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Figure 5: The scaling behaviour of (|-^ Y, x Tr ^o/^ll) f° r V < Vc (left), y > y c (right), 
at finite volumes and in a suitable range of The data is the same as in Fig. [|. The 
Ising model exponents are ~ 0.33 for y < y c and ~ —0.62 for y > y c (see, e.g., 
The continuous lines are fits to the 32 3 data. 



at y > y c is a well understood finite volume effect: for large enough volumes, it is 
proportional to 1/ V volume. 

At x = 0.30 we are substantially away from the first order region, and at y = y c there 
is a second order transition separating the two (Tr Afy ^ phases from the (Tr A^) = 
phase. Indeed, Tr A$ behaves exactly as magnetization in the 3d Ising model (this can 
be understood by comparing the regime x > Xtricriticai of the phase diagram in the left 
panel in Fig. |2| with that of the Ising model). We verify the Ising type of the transition 
by examining the critical behaviour of J2x Tr Aq/ g^\) in the neighbourhood of y c . 
For y < y c and in the infinite volume limit it approaches zero with the critical exponent 
(3 w 0.33: 

(\jpY, TlA o/9l\) oc [y c -yf. (29) 

When y > y c , (|i J2 X Tr ^o/flfl) = in the infinite volume limit. However, at a 
finite volume, its behaviour can be approximated in some range of y (where the linear 
extension of the lattice is still much larger than the longest correlation length) as 

<I^E T ^oA?3 3 D <xN-V\y-y c r</\ (30) 

where for the 3d Ising model 7 ~ 1.24. As can be seen from Fig. ||, the data agree 
rather well with these Ising model expectations. 

The measured (x, y)-plane phase diagram is shown in Fig. |6|. The first order line is 
converted into a second order transition at a tricritical point, which is approximately 
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Figure 6: The phase diagram of the 3d SU(3) + adjoint Higgs theory. The open 
symbols are results from the simulations, and the filled circle is the perturbative result 
in Eq. ([H]). The transition line is a polynomial fit to the data. 

located at x c « 0.26, y c « 0.104 (x Ci improved ~ 0.24 according to Eq. (g)). The fact that 
there is a second order transition at x > x c , is based on the observations that (1) the 
transition is not of the first order since the local order parameters behave continuously 
(Fig. (2) there is a real symmetry A — > — A which gets broken, and shows the 
scaling expected (Fig. |^), (3) the correlation length related to A seems to diverge at 
the transition point (Sec. |5.2| ). 

The critical exponents associated with the tricritical point in 3d assume their mean 
field values. We did not attempt to numerically analyze the critical properties at the 
tricritical point, since a meaningful analysis would have required simulations with much 
larger volumes and higher statistics than those used in this work. 

5.2 Correlation lengths 

The spatial correlation lengths, or inverse screening masses, were measured with a 
method using recursive levels of smearing and blocking of the gauge and Aq variables. 
We used up to 4 blocking levels, which gives 2 4 as the largest spatial extent of the 
operators. For technical details, we refer to [[|. 

In Figs. [7] and || we show the screening masses of the scalar operators Tr Aq and 
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Figure 7: The screening masses of the scalar operators Tr^Q and Tr A$ at x — 0.2 and 
0.3. In the broken phase, both operators couple to the same states. 
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Figure 8: The screening masses squared of the vector operators hi = Tr AqFu and 
hi = Tr AqFi2 at x — 0.2 and 0.3. In the broken phase, both operators couple to the 
same states. 

Tr Aq and the vector operators 

hi,i = e ijk Tr A F jk and h 2>i = e ijfc Tr A%F jk (31) 

at x = 0.20 and 0.30. In the symmetric phase, TrAg, Tr^Q have different quantum 
numbers and thus couple to different states; the same is true for hi, h^. However, in the 
broken A phase the operators can couple to each other, and should thus project to the 
same states. This is indeed observed in Fig. |7| within the statistical errors for x = 0.20. 
For x = 0.30 the signal of the Tr Aq correlator is quite noisy (and the errorbars do not 
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contain any estimates of the systematic effects), but the pattern should be the same. 
Note that this does not as such imply that there would be a discontinuity in the mass 
spectrum at y = y c , x = 0.30: it is just that in the broken phase Tr Tr A% couple to 
the same excitations, and to determine the mass of the first exited state would require 
a mixing analysis, such as in f2l |. In the symmetric phase, in contrast, TrA^Tr^Q 



automatically couple to different states. 

The scalar operator Tr Aq becomes "critical" at x — 0.30: it approaches zero at 
y = y c , as much as allowed by the finite volume. This is in accordance with the 
interpretation that Tr A^ acts like the magnetization of the 3d Ising model. A precise 
finite size scaling analysis of the critical correlation length is beyond the scope of this 
paper, and we rather interpret the right panel of Fig. ^ as a consistency check. 

In perturbation theory, there always remains an unbroken U(l) gauge symmetry 
in the broken phase. The vector operators hi and h 2 couple to the massless U(l) 
"photons" 7. However, these photons become massive due to the interactions with 
Polyakov monopoles. Since we expect the screening mass to be very light, we measure 
it at non-zero transverse momentum: we use the operator 

Oi (*s) = E T±A Fue a *'i' N (32) 

Xl,X2 

and correspondingly for h 2 (on the lattice, Try4o-Fi2 — > Tr A0U12). The screening mass 
is then measured from the asymptotic behaviour of the correlation function, 



(Oi(x 3 )O 1 (0)} oc exp(-a;3 V /(2vr/iV) 2 + m2 i ) . (33) 

The results are shown in Fig. [5[ When x = 0.2, the transition is strong and, within the 
statistical accuracy, the masses fall to zero in the broken phase. However, at x — 0.3 
the masses remain finite even in the broken phase. Due to the analytic continuation 
this implies that also at x — 0.2 the masses are still finite; they are just too small to 
be observed. 

It is interesting to note that the vector masses seem to vary rather smoothly when 
going from the symmetric to the broken phase, despite the fact that there is a true 
symmetry breaking phase transition. In fact, their behaviour is quite similar to that 
in the SU(2) + adjoint Higgs theory [|J, where no transition at all is observed at 
large x 



6 Conclusions 

We have numerically determined the phase diagram of the 3d SU(3) + adjoint Higgs 
theory. We find a symmetric phase with (Tr Aq) ~ 0, (Tr Aq) = 0, and a broken phase 
with (Tr Aq), (Tr Aq) non- vanishing and large. These are separated by a transition line 
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which contains a first order regime, a second order regime, and a tricritical point in 
between. We did not observe any other types of phases. 

From the statistical physics point of view, these results are of interest as a new quali- 
tative class in the types of critical behaviours that have been found in 3d SU(iV)+Higgs 
theories. (It should be noted that even more classes could appear when external 
constraints such as a magnetic field are added, or when the Higgses are replaced by 
fermions [^5] - the latter case is not relevant for finite temperature 4d theories, though. 
It can also be noted that a similar tricritical structure as found here can arise when 
there are several (fundamental) Higgses in the theory, related by a global symme- 
try [p6| .) It would be interesting to apply methods similar to those in [11 to a more 
precise study of the properties of the tricritical point. To our knowledge, the universal 
forms of the different probability distributions at the tricritical point have not been 
previously determined numerically in three-dimensional theories. 

From the QCD point of view, we re-emphasize the fact that even though the SU(3) + 
adjoint Higgs theory has the phase structure given in this paper, only the symmetric 
phase is an accurate effective theory for the 4d finite temperature theory, permitting, 
say, the determination of static correlation functions. In particular, the phase transition 
of the effective theory is not that of the 4d theory. (Nevertheless, it is amusing to note 
that the location of the tricritical point x c corresponds to T pa 0.57Aj^g according to 
2-loop dimensional reduction, not far from the 4d transition temperature T c pa 1.03Aj^g 



27].) One may also note that in the context of finite T QCD it has been suggested |2l 
that the critical line in the (T, /i)-plane for two massless flavours would have a similar 
tricritical structure. This clearly is not related to the tricritical structure studied here. 

Finally, from the point of view of the finite temperature SU(5) theory, the present 
SU(3) case shows that only some of the possible broken phases are dynamically realized. 
It should be interesting to see how this statement is modified in the SU(5) case, where 
there is a wider spectrum of possibilities. 
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